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Abstract 

We study the production of cosmological axions in the standard scenario in which a 
global string network forms at the Peccei-Quinn phase transition. We make detailed 
calculations of the axions produced by string loops, comparing these with estimates of 
other contributions from long strings and domain walls. We delineate key uncertainties 
in string network evolution, the chief of which is the magnitude of radiative backreac- 
tion. We conclude that axions produced by these topological contributions provide the 
significant cosmological constraint on the symmetry breaking scale / a and the axion 
mass m a . 

1 Introduction 

The nature of the dark matter of the universe remains an outstanding question facing 
cosmology today. Of the possible candidates, the axion is one of the most promising. In 
axion models an extra global £7(1)pq symmetry is added to the standard particle physics 
model to solve the strong CP problem [1]. This symmetry is broken at a high energy 
scale / a and the resulting pseudo-Goldstone boson — the axion — subsequently acquires 
an extremely small mass through 'soft' instanton effects at the QCD phase transition [2]. 
Because of substantial redshifting before it becomes massive, the axion is an ideal cold 
dark matter candidate. Viability, however, is dependent on compatibility between the 
predicted cosmological axion density and observational and nucleosynthesis constraints 
on the total and baryonic densities [3]. These constraints provide a significant upper 
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bound on the symmetry breaking scale / a which is near conflict with astrophysical bounds 
from red giants and SN1987a, / a > 10 9 -10 10 GeV [4]. 

If the universe passes normally through the Peccei-Quinn phase transition, a network 
of global strings will form by the Kibble mechanism [5]. (We note here that we will not 
be considering the alternative string-less scenario in which this transition occurs before 
an inflationary epoch — see, for example, ref. [6].) The decay of this string network into 
axions provides the dominant contribution to the overall axion density fi a [7]. However, 
the dependence of fi a on the energy scale / a has been the subject of much debate because 
radically different string radiation mechanisms have been proposed [7,8,9]. In a recent 
publication [10], we investigated the string radiation spectrum in great detail, using 
both analytic and numerical techniques. These findings were in broad agreement with 
the original work of Davis et al. [7,9] and contrary to the predictions of Sikivie et al. [8]. 
Taken at face value these results support a string constraint which would rule out the 
axion in the standard scenario. However, we also noted that, although the underlying 
physics of the earlier work was correct, the model employed for string evolution was too 
simplistic. 

In this letter we present a model for the evolution of an axion string network based 
on a marriage of numerical results for local strings [11,12] and the now- validated analytic 
radiation calculations for global strings which employ the Kalb-Ramond action [13,14,10]. 
Using this model, we calculate the density of axions radiated by string loops from the 
scale-invariant network properties (a previously overlooked source [15,16]). Using es- 
timates for these network parameters, the loop contribution to the axion density is 
compared to previous estimates for the contribution from long strings, domain walls and 
the homogeneous zero momentum mode. 

2 The evolution of an axion string network 

The evolution of axion strings is qualitatively very similar to the evolution of local 
strings due to their dynamical correspondence — as demonstrated numerically [10]. The 
additional long-range Goldstone field acts primarily to renormalize the string tension 
and energy density, 

/i^27r/ a 2 ln(t/«5), 
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where the string core width is 8 ~ J" 1 and we assume the typical curvature radius of 
the strings at a time t is R ~ t. Quantitatively on small scales £«t, global strings are 
affected by enhanced radiation backreaction; typically in a cosmological context axion 
radiation will be three orders of magnitude stronger than the weak gravitational radiation 
produced by local strings. This difference will alter small-scale features such as string 
wiggliness and loop creation sizes, but not the more robust large-scale network properties. 
As observed in numerical simulations of local strings [11,12], large-scale properties are 
remarkably independent of small-scale effects. As a first approximation, therefore, we 
employ parameter values found by these simulations. 

After formation, global strings experience a significant damping force due to the 
relatively high radiation background density. This frictional force eventually becomes 
negligible and the strings begin to oscillate relativistically and to radiate axions from the 
time U ~ 10- 20 (/ a /10 12 GeV)- 4 sec. 

From onwards, we can expect the axion strings to approach a scale-invariant 
regime in which the network properties remain constant relative to the horizon size 
djj = 2t. The overall density of strings splits neatly into two distinct parts, the long 
strings with length £ > t and small closed loops with £ < t. The long strings have an 
overall energy density 

Poo = C^, (1) 

where the correlation length scale of the Brownian network is £ (~ 1 ^ 2 t where ( 13. 

To describe scale-invariant loop creation and decay we must define several further 
parameters: First, we employ a to characterize the average loop creation size, that is, 
(£) = at. Secondly, the loop backreaction parameter n describes the radiation power 
from loops which is given by [13] 

P = YJl = K p , (2) 

where T a is a factor dependent on the loop trajectory, but not its size, which is estimated 
to be (r a ) ~ 65 [13,17] (exploiting similarities with gravitational radiation). As the loop 
decays into axions, its radius shrinks linearly 

e = e i - K ^t- ti) , (3) 
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where t\ is the loop creation time and i\ = l{t{). Typically, in a cosmological context 
we have k ~ (r a /27r)[ln(t/<5)] _1 0.15. Finally, we define the long string backreaction 
parameter 7; lengthscales below are smoothed by radiative damping in one Hubble 
time. Naively, one might expect 7 k, but the study in ref. [10] indicated that long 
string radiation, primarily in the second fundamental mode, was somewhat weaker with 
7 ~ 0.1k. The significance of 7 is that it should set the minimum loop creation size, that 
is, we expect 7 < a < n. Given these assumptions and energy conservation considerations, 
our scale-invariant model implies that the number density of loops in the interval i to 
i + d£ is given by 

vdi 

-('•«"« = whftw^ ' e ~ at ' (4) 

where v 0.40(a 1//2 . 

Near the QCD phase transition the axion acquires a mass and network evolution 
alters dramatically because domain walls form [18], with each string becoming attached 
to a wall [5]. Initially, the mass is temperature-dependent, but it only becomes significant 
when the Compton wavelength falls inside the horizon, that is, m(t)t ~ 0.75 at the time 

t ~ 8.8 x 10~ 7 A 2 ( -A V' 36 ( %r~) ^ ( — V ^ sec > (5) 

Vio i2 Gevy v 6xl0 " eV y V 60 J 

where A is a constant of order unity which quantifies parameter uncertainties at the 
QCD phase transition*, 

/ - \ °- 82 / A \ ~ - 65 / \f \ -°- 42 

A = ( m a ) ( A QCD \ ( A/qcd \ 



(6) 

6 x 10" 6 eV ) V 200MeV ) V 60 ) ' 
with the final axion mass m a = m a (/ a /10 12 GeV) _1 and the mass at temperature T given 
by m(T) = 0.1m a (A QCD /T) 3 - 7 [19]. 

Large field variations due to the strings collapse into domain walls at t. Subsequently, 
these domain walls begin to dominate over the string dynamics when the force from the 
surface tension becomes comparable to the tensional force due to the typical string 
curvature a ~ /i/t, 

/ r \ 0.36 / _ \ -2 / ,r \ 0.5 

-6 A 2/ Ja \ / \ /A/qcd 



t w ~ 1.7 x 10" b A 2 — ^ ^— '-^^ sec . (7) 

1 10 12 GeV/ V6xlO- 6 eV/ V 60 ' 



* Note that to calculate t, we assume an effective number of massless degrees of freedom Af in an epoch 
when its actual value is falling. 
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The demise of the hybrid string-wall network proceeds rapidly [5], as demonstrated 
numerically [15,20]. The strings frequently intersect and intercommute with the walls, 
effectively 'slicing up' the network into small oscillating walls bounded by string loops. 
Multiple self-intersections will reduce these pieces in size until the strings dominate the 
dynamics again and decay continues through axion emission. 

3 Calculation of the loop contribution to the axion density 

Given the loop distribution (4), we can calculate the energy density of emitted axions. 
The radiation spectrum from a loop of length £, averaged over various loop configurations, 
is given by 

= f 2 Jg(^), (8) 



dto 

where the function g{x) is normalised by 



oo 



g(x)dx = T aj , (9) 

'0 

and T a is defined in (2) (this approximates the loop spectrum which is actually discrete). 
We shall assume that loops are at rest, because any initial velocity will be redshifted and 
the net error when averaged isotropically over all loops should be relatively small. 

The energy density of massless axions emitted at time t\ in an interval dt\ with 
frequencies from u\ to u\ + du\ is 

r-OO 

dp^h) = dtxduxfl d£n{t,t 1 )tg{ew). (10) 
Jo 

Assuming Af constant, the spectral density can be calculated by integrating over t\ < t, 
taking into account the redshifting of both the frequency, u = a{t\ )/a(t)cu 1 , and the 
energy density, /? a oc a -4 . Neglecting the slow logarithmic dependence of the backreaction 
scale we have 

^ (4)= ^jC*'(T)^r*(r^ s [ ( ' /4,)i/ju< ] • (n> 

Under the substitution x = £/ti , z = ux(tti ) 1//2 , the range of integration is transformed 
and (11) becomes [16] 



2,, faut 



^ (t) = w7^i dzg{z) 



1 Z ^ _3//2 i\ a ^ _3//2 

UJKt/J V K 



(12) 



since the contribution from the lower limit can be shown to be zero for the range of u 
under consideration. This implies that the peak contribution to the axion density comes 
from those axions emitted just before wall domination. 

One can approximate the integrals of g{z) by noting that the dominant contribution 
comes in the range 4tt < z < 47m*, where n* is the mode beyond which the radiation 
spectrum of loops can be truncated due to backreaction. Assuming 47m* << unt and 
using the normalisation condition (9), the integral (12) becomes 

a\- 3 / 2 " 



1- (1 + 

K 



(13) 



did 3iOK 3 / 2 t 2 

From this expression we can obtain the spectral number density of axions dn^/dco = 
dp^l ' du . Integrating and comparing with the entropy density of the universe, s = 
27r 2 A/"T 3 /45, the ratio of the axion number density to the entropy at t w can be calculated 

as 



n 



6.7 x 10 c 



s 



a 

1- (1 + - 

K 



-3/2 



\ ji m a \ ( fa. 
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6 x 10- 6 eVy V 1012GeV , 
using typical parameter values T a 65, v 0.40Ca 1//2 and ( 13, Assuming number 
conservation after t w and using the entropy density so = 2809(To/2.7K) 3 cm -3 and crit- 
ical density p cr i t = 1.88 X 10 _29 /i 2 gcm -3 at the present day, one can deduce that the 
axion loop contribution is 



10.7 

. K 



a\-3/2 

1- (1 + - 

K 



where the Hubble's constant at the present day is H = 100/ikms 1 Mpc \ 0.35 < h < 
1.0, 

4 Other contributions to the axion density 

The contribution from long strings was roughly estimated in ref. [10]. The basis for 
this calculation was the radiation power per unit length for a typical sinusoidal per- 
turbation, dP/d£ 7r 3 / 2 /167t*, with the long string backreaction scale given by 7 ~ 
(7r 2 /8)[ln(t/(5)] _1 . Assuming the radiative dominance of this smallest scale (as ob- 
served in ref. [17]), one can calculate the spectral density of axions from long strings 

dp. vr 3 / a 2 ( 



dco 8jcot 



2 ' 



(15) 



This also assumes that the energy lost by the long strings doesnot interfere with the scaling solution 



Using similar methods to those used for loops 

n...*1.2fc-'A(^) , ( I5 A_) M ', (16) 

which, as before, is found to be roughly independent of the actual backreaction scale 
7. The considerable uncertainty of (16) must be emphasised given its sensitivity to the 
amplitude of small-scale structure and the overall long string radiation spectrum. 
A comparison of the two contributions (14) and (16) yields, 

a\- 3 / 2 " 



/a\3/2 

' ^8.9 



1- (1 + 

K 



(17) 



Either of the contributions could be dominant for the expected parameter range, 0.1 < 
a/n < 1 with equality at a/ 0.45. 

An order-of-magnitude estimate of the demise of the string/domain wall network [21] 
indicates that there is an additional contribution O a;( j w ~ (^ann/^) 3//2 (/a/10 12 GeV), 
where t 

ann ~ is the time of wall annihilation. This 'domain wall' contribution is 
ultimately due to loops which are created at the time ~ t w . Although the resulting loop 
density will be similar to (4), there is not the same accumulation from early times, so it 
is likely to be subdominant relative to (14). Both the long string and domain wall con- 
tributions will serve to strengthen the loop bound (14) on the axion; they are currently 
being studied in more detail [22]. 

5 Discussion and conclusions 

The estimates for the loop and long string contributions to the axion density are sum- 
marized in fig. 1 as a function of the relative loop creation size a/ n. This ratio expresses 
the key uncertainty arising from our inadequate understanding of long string radiative 
backreaction, 7 If we take the value implicitly assumed by most previous au- 

thors, a/fi w 1, then requiring fi a < 1 in the loop bound (14) we obtain a stringent 
constraint on the symmetry breaking scale, 

/ a < 6.0 x 10 10 GeV m a > 100/xeV, h = 0.5 , (18) 

or / a < 1.9 X 10 n GeV,m a > 31/xeV, for h = 1.0 (we have not included the parameter 
uncertainties of (6)). In this case, the axion is left a very narrow window which may be 
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Figure 1: The loop bound (14) and the (constant) long string bound (16) on the symmetry 
breaking scale f & for various values of OlJ K. The solid line corresponds to h = 0.5, the dotted 
line to h = 0.75 and the dashed line to h = 1.0. 

closed by the long string and domain wall contributions. However, as our recent investi- 
gations indicate [10], a parameter value as low as «//t ~ 0.1 may be more appropriate, 
that is, a regime in which the long string bound (16) is more important. This appears 
to provide an overall (uncertain) upper limit on the symmetry breaking scale, 

/ a < 2.6 x 10 11 GeV m a > 23/xeV, h = 0.5 , (19) 

or / a < 8.3 X 10 n GeV , m a > 7/xeV , for h = 1.0. Note that the weakest string bound is 
stronger than the early homogeneous zero-momentum estimates, / a <10 12 GeV [3]. 

We have summarised a model for string evolution and radiative backreaction which 
should clarify and correct the methods by which axion string constraints are calculated. 
Considerable uncertainties remain and it is clearly a matter of some priority to resolve 
the axion string radiative backreaction issue [22]. Nevertheless, we conclude that ax- 
ion emission by strings provides the key cosmological constraint on / a and m a in the 
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standard scenario. A window remains for experimental searches for the axion above the 
astrophysical bound / a >10 9 -10 10 GeV, but it is tightly constrained. 
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